Strong coupling of light with one- dimensional quantum dot chain: from 

Rabi oscillations to Rabi waves 
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Interaction of traveling wave of classic light with ID-chain of coupled quantum dots (QDs) in 
strong coupling regime has been theoretically considered. The effect of space propagation of Rabi 
oscillations in the form of traveling waves and wave packets has been predicted. Physical interpre- 
tation of the effect has been given, principles of its experimental observation are discussed. 
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Introduction. - Rabi oscillations are periodical transi- 
tions of a two-state quantum system between its station- 
ary states in the presence of an oscillatory driving field, 
see e.g. [T]. First observed on nuclear spins in radio- 
frequency magnetic field (3], the Rabi oscillations then 
were discovered in many other two- level systems, such 
as atoms exposed to electromagnetic wave 0|, semicon- 
ductor QDs [1], Josephson qubits jH], spin-qubits [^, and 
between ground and Rydberg atomic states [3|. Besides 
the fundamental interest, the effect of Rabi oscillations is 
promising for realization of binary logic and optical con- 
trol in quantum informatics and quantum computing. 

Complication of physical systems where Rabi effect is 
observed imposes additional features on the ideal picture 
[H of this effect. They are the time-domain modulation 
of the field-matter coupling constant (1, 0|, the phonon- 
induced dephasing (lol . and the local-field effect (l^ . 
Tsl . 14 1 - just to mention a few. New capabilities app ear 
in systems of two coupled Rabi oscillators 15, 13, 12, isl . 

30, El. 

In spatially extensive samples comprising a great num- 
ber of oscillators, the mechanism giving rise to Rabi os- 
cillations induces also a set of nonstationary coherent op- 
tical effects, such as optical nutation, photon echo, self- 
induced transparency, etc. iH]. This is because the sam- 
ple size exceeds significantly wavelength and propagation 
effects come into play. In low-dimensional systems prop- 
agation effects are also manifested but their character 
changes qualitatively. For example, the computational 
model of the coherent intersubband Rabi oscillations in a 
sample comprising 80 AlGaAs/GaAs quantum wells ji^l 
predicts the population dynamics to be dependent on the 
quantum well position in the series. This result demon- 
strates strong radiative coupling between wells and, more 
generally, significant difference in the Rabi effect picture 
for single and multiple oscillators. In the present Letter 
we build for the first time a theoretical model of a dis- 
tributed system of coupled Rabi oscillators and predict 
the new physical effect: the propagation of Rabi oscilla- 
tions in space in the form of traveling waves and wave 
packets. 

Model and equation of motion. - Consider an interac- 



tion of an one-particle excitation in an infinite periodical 
ID chain of identical coupled QDs with electromagnetic 
field. A p-th QD is considered as a two-level system with 
\bp) and \ap) as ground and excited states, correspond- 
ingly, and the transition frequency luq. Dephasing and 
dissipation processes inside the QD are further neglected. 
The coupling may originate from different physical pro- 
cesses (electron tunnehng, dipole-dipole interaction, etc.) 
and is accounted for in the tight-binding approximation, 
i.e., is assumed to be restricted to neighboring QDs. Let 
the QD chain be exposed to a plane wave traveling along 
the chain, E ~ e'x.p[—i{ujt — kx)]. The wavenumber sat- 
isfies the condition ka < 1 with a as the chain period, 
providing later on the continuous Hmit transition. 

In the one-particle basis, the Hamiltonian of the system 
can be represented hy H = Hq + AH, where the term 



Ho 



hujQ 



^<T+e*("''''^-"*) +H.C. (1) 



describes Rabi oscillations in non-interacting QDs. Here, 
(Tzn and are the FauH matrices for n-th QD, and is 
the real-valued Rabi frequency. The term AH accounts 
for the QD-coupling and has the form as follows: 



AH = ^ J2 iWn) (a„+p| + K) {bn+p\) + H.c, 

n p=±l 

(2) 

where ^ is the coupling constant. Equation of mo- 
tion has the form of one-particle Schrodinger equation 
ihdt l^*) = -ff 1^); the wave function can be written as 
l*(0> = Ep (Mi) K) + Bp{t) \bp)). Taking into ac- 
count (HJ) and |[2|), we reduce the Schrodinger equation 
to a system of differential equations, which directly cou- 
ples Ap with Bp,Ap±i and Bp with Ap,Bp±i. Carrying 
out then the continuous limit transition for the variable 
Ap{t) by Ap{t) A{x,t), Ap+i{t)+Ap-i{t)-2Ap{t) 
a'^dlA{x,t) and in the same manner for the variable 
Bp{t), in the rotating- wave approximation [l] we arrive 
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at the system of equations as follows: 

dtA = -l(a;o ^Be^('=--*) +z^a2a2A, (3) 

dtB = ^(^0 + + !^Ae-^('="--*) + i^a^dlB , (4) 

Eqs. ^ and (0]) describe light - QD chain coupHng in 
the framework of formulated model. Because we are in- 
terested in the strong coupling regime, the quantity Q.r 
can not be considered as a small parameter and further 
analysis of these equations is carried out without recourse 
to the perturbation theory. 

Traveling Rahi waves. - Let us consider elementary so- 
lution of the system l[3l)-|[4]) in the form of traveling wave: 

A ~ gi{h+k/2)x ^-i{u+uj /2)t ^ ^i{h-k/2)x^-i{v-uj/2)t 

where /i is a given wave number and v is the eigenfre- 
quency to be found. Solving characteristic equation of 
the system ©-([I]) with respect to v determines the eigen- 
frequencies of system by 



= Ca'/i' - * ± 



(5) 



where A tjQ — w is the frequency detuning, V = ^fca^ 
and <& — 2^(1 — k'^a? /?>). These two solutions correspond 
to two eigenmodes, given, respectively, by 

" " 2(1/1 + $ - A/2 - y/i - Ca2/i2) ' 
Bi{x,t) = c^e*(''-fc/2)xg--»(i^i-^c:j/2)t ^ 



and 



(J^^j^Qi{h-k/2)x ^-i(v2-uJ /2)t 



(7) 



where Ci_2 are normalizing constants. Either of these 
modes is a superposition of ground and excited states, 
whose partial amplitudes oscillate both in time and space. 
Binding of ground and excited states is caused by inter- 
action of light with QD chain and vanishes in the limit 
of rifl, 0. In that case, Eqs. ^ and lO describe 
QD-chain excitons in equilibrium and inverse states, re- 
spectively. Retaining in the expansion terms linear in 
\^r\ and simultaneously substituting A ^ A — iO we ar- 
rive at the intermediate case of excitons weakly coupled 
with electromagnetic field. 

Space oscillations of the partial amplitudes are due to 
QD-coupling and vanishes in the limit of ^ 0. Thus, 
each of these modes can be interpreted as a Rabi wave 
with the frequency determined by Eq. In general 

case these waves are excited simultaneously, while any 
of them can be excited separately by a proper choice of 
initial conditions. 




FIG. 1: Typical dispersion law for Rabi waves. ^ = 

A — O.SSlfl, ka — 1. Turndown of ha corresponds to the 

range of applicability of the continuous limit model. 



Typical dispersion characteristics of Rabi waves are de- 
picted in Fig. [TJ It should be noted that at given ampli- 
tude and frequency of external field their frequencies have 
continuous spectra (wave number h varies continuously) . 
Dispersion dependences depicted in the figure are asym- 
metric: i'i,2{h) 7^ i'i,2{—h). Physically, this is due to the 
presence of preferential direction, which is determined by 
the direction of light propagation along the chain. The 
QD coupling leads to the inequality vi ^ —V2. That 
is why, unlike to single QD, in QD chains the inversion 
oscillates anharmonically. These oscillations can be rep- 
resented as amplitude-modulated harmonic oscillations 
with the frequency {vi — 1/2) / 2 , while the modulation fre- 
quency is given hy {vi + V2) /'i. In traveling Rabi wave, 
the inversion is constant in space. Because vi^2{h) are 
real at any real h, the system is stable j23|. Now, let us 
analyse the dispersion characteristics hi ^2(1^), assuming 
the frequency to be a given parameter. It is seen from 
Figdithat in some frequency range hi is complex and 
is real for real ly. It corresponds to non-transmission of 
the first Rabi wave [23|. The frequency range in which 
both of /ii,2 are complex also exists. This case corre- 
sponds to complete non-transmission of the Rabi waves 
with given frequency. 

Note that the eigenmodes ^ and ^ each comprise 
traveling waves with different wave numbers h ± k/2. 
Physically, this means that the Rabi wave propagates in 
an effective periodically inhomogeneous medium formed 
by spatially oscillating (with period 27r/fc) electric field. 
Therefore, the diffraction is developed in the system. In 
the limit /c — > the medium turns homogeneous and the 
diffraction effect vanishes. 

Assuming the frequency to be a given parameter and 
solving Eq. ^ with respect to wave number, we obtain 
for fc = 0: 



hi,2 = ± 



1 



1/2 



(8) 



where external signs ± correspond to two directions of 
propagation while signs ± before internal radical cor- 
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respond to two types of Rabi waves indexed by 1 and 
2. It directly follows from ([H that electric field forms 
an effective medium for propagating Rabi waves. In in- 
homogeneous electric field the Rabi frequency becomes 
coordinate-dependent, ^Ir = ^Ir^x), and therefore the 
medium becomes inhomogeneous too providing refiection 
of Rabi waves and their mutual transformations at the 
inhomogeneities. In that way one obtain a unique ability 
to control the processes of the refiection and dispersion of 
Rabi waves by varying the light spatial distribution. As 
a potential realization scheme we indicate the interaction 
of QD-chain with Gaussian light beam (or superposition 
of such beams) with the beams' widths and mutual dis- 
position as controllable factors. 

Rabi wave packets. - The process of the excitation 
transition opens up new opportunities for controlling the 
dynamics of Rabi oscillations. For identification of con- 
trol factors we need to know general solution of the 
system ll3l)-|[4l). To find it, we first introduce the vari- 
ables u{x, t) = A{x, t) exp [i(ujt - kx - 2$t)/2], v{x, t) = 
B{x,t)ex\)[—i{ujt~kx + 2'^t)/2]. For these variables, 
Eqs. ll3])-(l4]) are reduced to the form as follows: 

dtu + -I- VdxU — i^a^d^u — ^^^v — , 

Otv — VOxV — i^a o^v —u = . 

This system can be solved exactly by using the Fourier 
transform with respect to x. Finally we arrive at 



u{x,t) = J [u{h)^-{t) + i{h)Mt)] 6^"^"-^° '''^dh, (9) 

— OO 

oo 



where 



COST ± i{Ah/^h) SUIT, T 



A,, = A 



2Vh . The function u{h) is determined by the initial con- 
dition 



u{h) 



1 

2^ 



u{x, 0)e 



-ihx 



dx 



(11) 



and the function v{h) - analogously. It is easy to verify, 
that the solution obtained satisfies the probability con- 
servation law -I- \v{x,t)\^]dx — 1 for any 
t > Q. A typical space-time distribution of the spatial 
density of the inversion (the inversion per single QD) 
w{x,t) = a[|u(2:,i)p — |w(x,t)p] is shown in Fig. [2k- As 
follows from Eqs. (|9l)- (fTQl) the depicted space-time dy- 
namics of w{x,t) corresponds to the evolution of a Rabi 
wavepacket defined as a superposition of Rabi waves with 
continuous spectrum of Q.h. Physical interpretation of 



the picture predicted to observe in the QD chain can be 
given on base of the collapse- revivals concept [H . Distinc- 
tive feature of the case considered is that the distribution 
of collapses and revivals is permanently varied in space 
and time. 
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FIG. 2: Space-time distribution of the inversion in the 
QD chain (color online), a) A single Gaussian wavepacket 



A[x,Q) = exp(-x72cr2)/A 



B(a;,0) = 0, A = Vk. 



b) Two counterpropagating identical Gaussian wavepackets: 

A{x,0) = exp[-(x - 3ct)V2ct^]/a/4^, B(x,0) =exp[-(a; + 
3cr)V2(7^]/v^^, A = 0. 



In both cases ka = 0.33, a = 5a, 



Although variation of the inversion density, depicted 
in Fig. [21 in arbitrary point of the space occupied by 
the Rabi wavepacket is not too large, an integral char- 
acteristics presented in Fig. [3|- the "integral" inversion 
w{x,t)dx - of initially unpolarized QD-chain oscil- 
lates between -1 and 1, thus indicating presence of strong 
light-QD coupling. 

Note that oscillations of the integral inversion atV j^O 
damp with time (see Fig. [3]), whereas such a damping is 
absent at = and integral inversion oscillates har- 
monically in the range from -1 to 1 (dotted curve in Fig. 
[Sj. Such a behaviour indicates appearance of a specific 
mechanism of collective dephasing. Physically, this is 
because the effective detuning Ah and therefore the car- 
rier frequency flh of Rabi oscillations constituting the 
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wavepacket depends on h (Doppler shift). As different 
from that, the condition V = keeps the /i-dependence 
only in the ampUtude modulation frequency {i^i + z/2)/2 
and thus does not result in dephasing. In the weak cou- 
pling limit the indicated dephasing mechanism is analo- 
gous to the Landau damping in plasma. 

The dependence of the frequency of Rabi oscillations 
on V shows that the value A = is not optimal for the 
effect observation. Optimization of A allows increasing 
the intensity of Rabi wave and the dephasing time, see 
dashed curve in Fig. [3l In the frequency domain, fine 
tuning of the system to the resonance is achieved by the 
variation of V (changing the angle of incidence of Hght). 
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FIG. 3: Temporal dependence of the integral inversion at the 
ioput parameters as follows: A = 0, fc = (dotted line); 
A = 0, fea = 0.33 (solid line); A = Vk, ka = 0.33 (dashed 
line). In all cases, ^ = 3flR and a = 5a. 

Interaction of two counterpropagating identical Gaus- 
sian Rabi wavepackets elastically colliding at x = is 
shown in Fig. [2b- Although the inversion oscillates in 
time and moves in space, integral inversion of initially 
saturated QD chain (/T^ w{x, 0)dx = 0) does not expe- 
rience oscillations: this quantity equals zero for alH > 
and arbitrary values of flR. 

On experimental observability of Rabi waves. - The 
theory presented can be extended to a various physi- 
cal situations, such as quantum dynamics of an elec- 
tron in QD chain. Rabi waves are realized via optically 
induced transitions between the size-quantized electron 
levels. Transition of electron from one QD into another 
occurs by means of tunneling through the potential bar- 
rier 2l|, |25||. Another example is the Rabi oscillations of 
two-electron entangled state (|01)-|-|10))/ V2 taking place 
in two neighboring QDs due to dipole-dipole interaction 
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20J. The model developed describes 
the motion of this two-electron state as a single whole; in 
this case the wavefunction |^) is the envelope function. 

Theoretical analysis carried out has shown that the op- 
timal excitation of Rabi waves require the coupling fac- 
tors of both neighboring QDs and single QD with field 
to be comparable by magnitude: ^ ~ Another nec- 
essary condition being imposed on the Rabi frequency is 
essential exceeding over all intrinsic relaxation rates. For 



typical QD structures Q, this condition is satisfied for 
Rabi frequency varied over a wide range hflu 0.001 — 1 
meV. This range corresponds to reaHstic electric field 
variation E ~ 10^ — 10^ V/cm. As a consequence, the 
interdot coupling constant ^ also varies between 1 ^eV 
and 1 meV what is practically achievable 15, 13, 21 1 for 
typical interdot distances 4 — 20 nm. 

Experimentally, the Rabi waves can be detected in res- 
onant fiuorescence spectra of spatially extensive samples 
by the presence of new spectral lines, additional to the 
Mollow triplet, as well as by the Doppler shift and broad- 
ening of the triplet lines, etc. Of course, highly ordered 
chains of uniform QDs are required to exclude nonho- 
mogeneous broadening, which may hide the effect. Im- 
pressive progress in growing of perfect nanostructured 
successions achieved in last years (e.g., see ji^) is very 
promising for that aim. 

Rabi waves can be observed in systems of another 
physical nature such as quantum electrical circuits 
if one proceed from two coupled Josephson cubits in mi- 
crostrip resonator |3] to a distributed structure of such 
elements imposed to interqubit interaction. In particu- 
lar, in that structure the Rabi wave frequency goes down 
to microwaves and the field intensity necessary for Rabi 
waves excitation decreases. 

Conclusion. - In this Letter, we have predicted the 
exitance of Rabi waves - wave propagation of the inver- 
sion in spatially extensive systems of coupled oscillators. 
The system has been exemplified by an ID-chain of cou- 
pled QDs exposed to an intensive traveling light wave. 
Spatial propagation of Rabi oscillations in the form of 
traveling waves and wave packets is shown to occur in 
the chain. The propagation is predicted to be accompa- 
nied by the damping of Rabi wave in time manifesting 
new mechanism of collective dephasing, which is an ana- 
log of Landau damping of exciton-polaritons extended to 
the strong light-matter coupling. 
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